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ABSTRACT 


The stability of systems is investigated which can be 
described by differential equations consisting of a linear time- 


varying part and one nonlinear term. 


First, only the linear part is considered. Definition 
and properties of the impulse response are discussed in detail. 
Closed form solutions are evaluated for some special systems. For 
the general case, approximations are determined by comparing two 


integral equations. 


Then, utilizing the results of the linear part, an 
integral equation is set up for the complete system, including 
the nonlinear term. Sufficient stability conditions are derived. 
It is shown that the appropriate Volterra series is convergent under 


relatively weak restrictions. 
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CHAPTER I 


INTRODUCTION 


1.1 Preliminary Remarks 


Continuous systems with lumped parameters in general can 
be described by ordinary time-varying differential equations. Since 
it is difficult or often impossible to solve these equations, the 
investigation of a physical system usually requires admissible assump- 
tions which simplify the mathematical description. Fig. 1.A illustrates 
a customary procedure for linearizing a system or to make it time- 


invariant (stationary). 


For a linear time-invariant system, a complete mathematical 
theory is available which has been successfully applied to numerous 
problems. However, in those cases where the results, obtained from the 
linearized time-invariant system, are insufficient or the system does 
not allow such a simplification, nonlinearities and (or) time-varying 


parameters have to be included. 


In spite of the huge amount of work which has been done in 
the field of nonlinear time-varying differential equations, the results 
do not form a complete theory yet. This is mainly due to the fact that 
solutions, in closed form, exist only for a very few special classes of 
systems. Consequently, for the general case, one has to resort to 
various approximations depending on the particular problem and the type 


of system under consideration. Moreover, it is often preferred not to 
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Figure 1.A 


determine the complete solution, but to investigate only those prop- 
erties which have to be known for a certain application, for example: 


zeros, extrema, monotonic or oscillatory behavior, envelopes, etc. 


From these introductory remarks it follows that the objective 
of this thesis has to be limited to the investigation of some properties 


of a class of systems, which will be outlined in the following section. 


ieee. Obj ectiverandeeutline of the Method 


In this thesis, the following class of systems will be dealt 


with: 
iota G) et eeng (x,t) * r( ty. Clee) 


L(p,t) is a nonsingular, time-varying linear operator defined as: 


-l1 
Tipe) p> + a1 ft) p” Eee SES 


CDS) 
eS a, (t) pt a(t) 


where: 
ee 
es 
ECL) is ean sarbitrary driving function, 
€ is a positive constant, and 


o(x,t) is a nonlinear function, specified in detail in 


Chapter IV. 


Equation (1.2.1) describes, for instance, the feedback systems 


imekigeel.b and Big. 1sC. 
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The tunction “Ct)) in €1.2.1) is defined as: 


r(t) = L(p,t) r,(t) 


Figure 1.B 
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Figure 1.C 
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System (1.2.1) becomes more familiar, when the special 
case L(p,t) = L(p) is considered, in other words, in the case of 


the stationary linear differential operator: 


L(p) =p’ ta y eat agi appt a, (1.2.3) 
where Qon-1 ? °°* 2 F601 > B50 TE constants. 
Then (1.2.1) becomes: 
Lop) x Choc ec é(x.t) & r(t). (Giles 


One of the properties of (1.2.4), which has been the 
object of extensive research, is the stability for various inputs 
or initial conditions. The two best known results are the Popov 
and circle criteria, which both use a frequency domain approach to 
derive sufficient stability condittons: Brockett [1], Barrett [2], 
Trott and Christensen [3], and others apply certain results of 
functional analysis to equation (1.2.4) in order to get sufficient 
stability criteria and to generate a series of successive approxima- 


tions. 


A prerequisite for applying these methods is the determin- 
ation of the inverse operator ti i(p) » so that the solution of (1.2.4) 


can be written as: 
Se iiPrety tiger eG eae. (1.2.5) 


or in form of the convolution integral: 
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where: 


y(t) is the solution of 
(py y(t) = 0 


subject to the given initial conditions, and h(t-t) is the so 
called impulse response. 
For the time-invariant operator (1.2.3), the functions y(t) and 


h(t-t) can be obtained-in closed forn. 


If it is desired to use a similar approach for the inves- 
tigation of the time-varying system (1.2.1), one could formally 


rewriteeGlecwleas: 
x(t) = L+(p,t) {r(t) - © 6(x,t)} G@.2.7) 
or in integral form: 
ic 
(Ces ay Gb) er | Gtr) meiner )e=.e.0( xr) | edr Gie2ec) 
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The major part of this thesis will deal with the problem 
of justifying (1.2.8) and finding closed form solutions and approxima- 


trons foOneey(t)eeands h(t, 1) 


In Chapter II, basic properties of linear time-varying 


differential equations are reviewed which are relevant to the impulse 


0 = (3)e (q)d 


md 


oe sda at (7+3)d Bas (ang tiboo> tle t2 gees . 
| 4 
bab (3)y amotsonud sz , (€.5.£) t03878q0 31 
\i0k beaols nb bentssdo . 3 
— << 
; a a 7 
“eovat oid “302 donoxags wlimta & seu o3° basteeb ab 3h. 23 os) F 
| eh tiaeo? Ebiew Sacks (1.8.1) medaye si eas 38 ett 3 iy 
othe | _ 


ne, 


2.) (806 9+ 3)a) (ag) e= Ge 
% d 7 : ‘ 


(8.5.1) « 7h [(r,x)o 3 - (alae ‘| my Gy = Gx 
9 oO. - 
| | a 


@oidorq od3 Widiw.ieeb ftw atesds eti to 32Bq - et 
eee ete net 


ie 


response h(t,t) . Several classes of systems, which permit solutions 
in closed form, illustrate the properties of h(t,t) and the manner 


in which it is computed. 


In Chapter III, a unified method is presented for approx- 


imating certain types of linear systems. 


In Chapter IV, the results of Chapter II and III are applied 
to (1.2.1). Sufficient stability conditions are derived and the conver- 


gence of the Volterra series, developed from (1.2.8), is proved. 
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CHAPTER II 


LINEAR TIME-VARYING DIFFERENTIAL EQUATIONS 


2.1 Introduction 


This chapter deals with the solution of linear time- 


varying differential equations. 


In section 2.2, basic properties of linear differential 
equations are reviewed, which are needed in order to derive h(t,rT) 
as the one-sided Green's function. In this manner, one arrives at 


a simple and rigorous definition for h(t,t) . 


In section 2,3, h(t,1) is evaluated in closed form for 
some special systems. A transformation in the time domain is intro- 
duced which reduces some time-varying systems to systems with constant 


coefficients. 


2.2 Review of Basic Properties [4] 


Consider the following linear differential equation: 


E(p,e) x(r)e=i (pea _,(t) file auras (2.2.1) 
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subject to the initial conditions: 


CA?) 


The solution of (2.2.1) is the sum of the homogeneous 


and the forced solution: 
rel Geal) CP AG cy ACen oe G2 23.) 
The homogeneous solution y(t) satisfies: 
Lipa) ~Co)e=20 (22254) 


andgincludessmie initialecondi tions (2.2.2), that is: 


oad) 


(n-1) > 
y (t)) * Oo 


‘ eee 6 Ga) = an ; (23255) 


Existence and uniqueness of y(t) are guaranteed by the 
Cauchy-Lipschitz condition which, in the case of (2.2.1), is 
satisfied if the coefficients a, (t) are continuous functions of t. 
Throughout this thesis it will be assumed that the coefficients a, (t) 
satisfy the more restrictive condition to have derivatives of any 


desirable order. 
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The homogeneous solution y(t) can be expressed as a 
linear combination of n linearly independent functions o, Ct) ; 


each of which satisfies (2.2.4). 


Thus: 
y(t) = c, o>, (t) dearly trtiaptin(e)ic. 


Linear independence implies that the constants c, can be determined 


from the initial conditiong+@ 32.5) 


The forced solution z(t) is defined by: 
mCp ec ez) ear C(t) Od) 
for the following initial conditions: 


need) x. peas (ee) ee z(t.) = 2(t.) 20. (2.2.8) 


The solution of (2.2.7) is conveniently written as convolution 


integral: 


i 
z(t) = | Ger ce iy ea el le Yebal y, (252.9) 
Ke 
O 


A simple and rigorous procedure to derive (2.2.9) and 
to define h(t,t) uses the method of ‘Variation of Parameters'. This 
way avoids in addition the lengthy introduction of a mathematically 
acceptable 6-function [5]. In order to illustrate this procedure in 


some detail, consider an example of second order: 
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z(t) + a(t) 2(t) + a(t) 2(t) = r(t) 
with the initial conditions: 


z(t.) = 0 and z(t.) = 0 


(242,10) 


(2s2rasy 


Let o,(t) and bo (t) be two linearly independent solutions 


satisfying the homogeneous solution: 
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Then the method of 'Variation of Parameters’ tries to determine 


two functions 6, (t) and Ey (t) such that the solution of (2.2.10) 


can be written as: 
a(t) = €,(t) $,(t) + E(t) $5(t) 
Differentiating (2.2.12) gives: 


The first equation chosen for the unknown functions an 


is: 
25, Ol Sd ama 
Consequently, the second derivative of z(t) is: 
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This yields the second equation for E, (0) and Ey (t) : 
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can be solved for © and ae » so that: 
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(22 AT) 


D(t) &, = + x(t) $4 (t) 
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Gece ot) 


Integrating (2.2.17) 
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and substituting (2.2.19) into (2.2.12) gives: 
z(t) = &,(t.) $,(t) + €,(t.) o,(t) 
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Observing the initial conditions (2.2.11) results in: 
6, (t,) = 0 and E(t.) = 0 
so that (2.2.20) becomes: 


ear (Tt) de (BAe Earby) 
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Cemparing) (2.2.21) with (2.2.9) yields the definition of. h(t,7): 
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for the initial conditions 


ht(tec) = 0 and p htt,t) =] (oe) 


An immediate consequence of (2.2.23) and (2.2.24) is that the exis- 
tence and uniqueness of h(t,t) is assured, whenever (Qe 2a) eas a 


unique solution which is quaranteed by the Cauchy-Lipschitz condition. 


Another property of h(t,t) which will be of interest 
later is the dependence on the parameter 1. For the second order 
example (2.2.10) it will be shown that if: h(t,t) is the impulse 


response of 
ate es vad a(t) y= 0 
and eg(t,t) is the impulse response of the adjoint equation: 
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a dt iL a fe) a 


which implies: 
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a | t= 
then 
h(t,t) = -g(1,t). (215-98) 


In order to prove (2.2.26), the following expression will be integra- 


ted by parts: 
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(lower case indices denoting partial derivatives) 
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or the desired result: 


= g{t,t) = h(t,7) 


If the previous derivation is extended to higher order 


systems, the following result is arrived at: 
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If t is the independent variable, h(t,t) satisfies 


the adjoint equation: (note: p oe ) 
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For the convolution integral (2.2.9), h(t,t) has to be determined 


on the intérvala.0 ace th< = & 


2.3 Solutions in Closed Form 


In this section, solutions in closed form will be obtained 
for some classes of linear homogeneous differential equations. By 
imposing the conditions (2.2.27), the appropriate impulse response 


will be obtained in each case. 
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Solution: 
t 
a a (ua) du 


y(t) = y(t ye ° 


According to (2.2.27), the impulse response is the 


solution of (2.3.1) which satisfies: 
y(t.) = y(t) = 1 


Consequently: 


-{a.(u)du 
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b. Second order systems 
y + a,(t) y+ a(t) y = 0 (e372) 


in order to ebtainethe solution of (2.3.2) in closed form, 
certain assumptions have to be made regarding the coeffi- 


cients a, (t) and a(t) 
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Then (2.3.2) can be rewritten as follows: 
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The general solution of (2.3.3) is: 
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For the special initial conditions (2.2.27): 
y(t.) = y(t) #0 and y(t.) = y(t) = 1 


the general solution (2.3.4) becomes the impulse response: 
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reduces (2.3.5) to a differential equation with constant 
coefficients. 


For the following specific numerical example: 


oe 4 é 9 
Veer ane yr (23.74) 
ttl Ce+1)° 


the substitution € = In(tt+l) yields the general solution: 


Ree 
2 I = nes (gt ne =< we 


dani: salognt sila siti meee > ee Pe 
F 3 
vb (ub (4) (6 | +e |» (7,3) pS 
Ps ' 


:ndigevps Istinsrstiib @'tslua rsbtened 4 


lion. 30 ae Gaver 
ese.) oe tae ee a 
gntwol lot oft .etnsdenos sitB 30" bas ais sti fe 
| : nora Siem = 
sae aes 
DES) (stra = 3 (3)m = (dy 


Nid - : a =). ' 


inetenod dyiw notseups lelinsratith 5 03 (2.6.9) eaadbes - 
‘i : a 
7 ? > Z 


ba 
4 


so | Pais age ues atlas a sits 


SBS) | O48 Xe signe 2. - 
eaihlenndits woh + 00 


eo aed 


cae ttl 
Yaa y Gti) Ally aelee ay 
(23357) 
J t +1 2 
eG 53) sarap aa Gus ae sp) 
from which the impulse response is obtained as: 
Th a2 
hn. Gtr) s 7) Cur at) (2 no) 


For time-invariant linear differential equations, the 
general solution and consequently h(t,t) can conveniently be 
evaluated by using the Laplace transformation. Time-invariance 


implies the following property: 


y(t) s subject to certain initial conditions at° t°= (0 


is equal to y(t - t.) for the same initial conditions 
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The solution of (2.3.9) can be obtained from: 
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The solution of (2.3.9) for the delayed initial conditions y(t.) ay 
and y(t.) ae is simply y(t - t) where y(t) is defined in 
(2.3.10). This implies for the impulse response: 
h(t,t) =h(t-t) , where h(t) = L7> [4=—+—— ] 
Since the solution of a time-invariant differential equation 
can easily be determined it would seem logical to search for a trans- 
formation which reduces a time-varying system to a system with constant 


coefficients. 


The substitution (2.3.7) for Euler's equation can be taken 
as one example of this kind of transformation. Muraview [6] proposes 
a more general method for a class of second order systems which will 


be dealt with in the remainder of this section. 
Suppose the following differential equation 
ad Se AES) y + a(t) y = 0 (53711) 


is defined for t > 1, and is to be solved subject to the initial 
conditions y(t.) and y(t.) . The main idea of this method is to 


introduce the following transformation: 
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n" +n’ [a,(t) a(t) - (t)] + nlaj(t) + aM(t)] = 0 (2.3.14) 
In order to reduce (2.3.14) to the following time-invariant system: 
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two constants m4 and m have to be determined such that: 
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Then (2.3.14) can be solved easily and the original solution y(t) 
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The method will be illustrated by the following example: 
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The constants m and m, can, for instance, be chosen as: 


Meeees, m, = 2 (223417) 


so that 


a(t) ve and Ee J] - - (2534-5) 


For the numerical values in (2.3.17), equation (2.3.15) has the 


general solution: 
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which, observing (2.3.12) and (2.3.18), can be rewritten in terms of 
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According to definition (2.2.27) the appropriate impulse response is: 
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As a concluding remark, the work by Kamke [7]*should be 
mentioned which contains an extensive list of differential equations. 
This reference might be helpful in finding a closed form solution for 


a particular differential equation. 
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CHAPTER III 


APPROXIMATIONS 


3.1 Introduction 


In Chapter II some classes of linear systems were discussed 
which could be solved in closed form. Unfortunately, these systems 
represent exceptional cases and should not conceal the fact that arbi- 
trary linear systems of higher than first order cannot be solved explic- 
itly. Therefore h(t,t) and consequently the convolution integral 
(1.2.8) cannot be obtained in closed form so that approximations have 


to be introduced. 


For this purpose,(2.2.4) will be transformed into an 
equivalent integral equation. This will be discussed in section 3.2. 
It will be shown that the sequence of successive approximations con- 
verges for the assumptions made in Chapter II. In section 3.3, a 
bound on the absolute value of the impulse response will be derived. 
Finally, in section 3.4, several examples will illustrate the methods 
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3.2 Linear Differential Equations and Linear Integral Equations 


A linear differential equation, subject to certain initial 
conditions, can as well be represented by a linear Volterra integral 
equation of the second kind. It should be pointed out that represent- 


ing the system by an integral equation does not imply a closed form 
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solution, and does not even imply an easier way to obtain the solution. 
However, as far as approximations are concerned, an integral equation 
offers a much wider range of methods than the original differential 


equation. 


In order to illustrate the transformation from a differen- 
tial equation into an integral equations in some detail, the following 
third order system will be used throughout this and the next section. 


Consider: 
y+ a, (t) y + a, (t) y + a(t) y = 0 (222015 
defined for t > ty and subject to the initial conditions: 
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In order to transform (3.2.1) into an integral equation, 


the following differential equation is introduced: 
yt Gace) y+ aye) y + ei(t) y = 0 (3798) 


where for the time being it will be assumed that the coefficients 
c(t) : c,(t), and Co (t) are "suitably" chosen and that (3.2.3) has 


a closed form solution. 


Adding (3.2.3) to both sides of (3.2.1) and transposing 


the terms with a, (t) yields: 
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Since the LHS of (3.2.6) was assumed to have a closed form solution, 


equation (3.2.6) can be expressed by means of convolution as: 
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and g(t,u) is the impulse response of (3.2.3) satisfying: 
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From (2.2.28) it follows in addition: 
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Integrating (3.2.7) by parts gives: 
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According to definition (2.2.27) the general solution 
(3.2.11) becomes the impulse response if the following initial condi- 


tions are substituted at ty pete a 


y(t) = 0, y(t) = 0, y(r) = 1 (ono) 


which gives: 
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h(t,t) = g(t,7) + | k(t,u) h(u,t) du (33S) 
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Equation (3.2.15) can be interpreted in the following way: 


The impulse response of (3.2.1) is equal to the 
impulse response eg(t,t) of the LHS of (3.2.6) 
plus a "correction" term (fs see du , due to the 


RHS of (3.2.6). 


The integral equations (3.2.11) and (3.2.15) are equivalent 


to the differential equation (3.2.1) for the initial conditions (3.2.2) 
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and (3.2.14) respectively. As mentioned earlier, the integral repre- 
sentation does not imply an easier way to a possible explicit solution. 
However, it is a prerequisite to prove convergence of successive iter- 


ations if the solution of (3.2.1) is to be approximated in this manner. 


The standard approximation method, sometimes called 


"Picard method of successive approximations", proceeds as follows: 
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increases, it must be proved that: 
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Chu and Diaz [8] use an elegant method to prove the conver- 
gence of (3.2.20) which will be reproduced here in a slightly differ- 


ent way. 
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Equation (3.2.34) can be interpreted as follows: The 
series (3.2.34) converges uniformly on every finite interval [0,T] , 
provided (3.2.26) is satisfied. This implies that y(t) can be 


approximated by a finite number of terms s(t) such that the error 


New 
ly(t) - )k s(t) | does not exceed a given ce . Therefore, from 
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the convergence point of view, the iterative approximations, defined 
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However, if the terms s,(t) are to be determined from 
(3.2.16) and (3.2.18), one will find in many cases that the evaluation 
generates too voluminous expressions or integrals which cannot be com- 
puted in closed form. Due to this difficulty, the result of the last 
proof can in general not be exploited to obtain an arbitrarily accur- 


ate approximation. 


S25 Integral Inequality 


A less demanding, but for later purposes still acceptable, 
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yields the following inequality: 
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While (3.2.11) defines a unique function, the inequality (3.3.1) is 


valid for a set of functions, one of which is ly(t)| . 
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the equality sign holds in (3.3.1) and the bound on ly(t)| is 
identical to y(t) . This is of course the ideal case. Therefore it 
is desirable to transform the given differential equation into an 
integral equation which has the properties specified in (3.3.2). One 
possible way of achieving this may be accomplished by the choice of 
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aspect of a suitable choice can be specified as follows: "Choose the 
coefficients so that the resulting integral equation is not changed 


"too much’ by taking the modulus". 
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the following property is true: 


y(t = ¥tt) for 0:2 t:< T (35955) 


In order to prove this seemingly self-evident fact, some basic 


properties of linear integral equations are reviewed. 
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Utilizing (b) and (c), Theorem I is now readily proved. 
Because of (c) every ly(t) | which satisfies (3.3.3) is bounded 


from above by: 
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In order to apply Theorem I, Y(t) has to be determined 
in closed form. The simplest form of K(t,u) which allows a closed 


form solution is: 
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so that (3.3.4) can be simplified to: 
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Equation (3.3.10), together with property (c) represents 
a generalization of various forms of the Bellman-Gronwall lemma which 
can be found in Hsu and Meier [10], Sansone and Conti [11], and 
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Inequality (3.3.12) can easily be verified by the methods 


developed here. 
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The appropriate integral equation (3.3.4) can conveniently be solved 
by using Laplace transform (see Doetsch [13]). Taking the Laplace 


transformation of: 
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Two special cases which will be of importance for the next 


section will be now discussed in more detail. Consider: 
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Differentiating (3.3.16) mn times gives an n-th order differ- 
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Equation (3.3.20) is a differential equation with constant 
coefficients. Therefore Z(t) can be solved in closed form. 
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Before the methods, discussed in the last two sections, 
are applied to certain examples, the basic procedure will be reviewed 


briefly. 
The differential equation (3.2.1) is given for y(t) . 


Equation (3.2.1) is converted into a suitable integral equation 
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Inequality (3.3.1) is compared with the integral equation (3.3.4) 


which satisfies the conditions of Theorem I. 


The comparison gives the desired result: |y(t)| se YOUR Ys 
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In this section, several examples will illustrate the methods, 
which were described in general in section 3.2 and 3.3. The efficiency 
of these methods depends to a large extent on equation (3.2.3) which 
has to satisfy two requirements. Firstiy, it must have a closed form 
solution, secondly, the corresponding integral equation must yield an 
acceptable bound on the appropriate solution. Unfortunately, there is 
no general rule according to which (3.2.3) could be determined so that 


the previous requirements are satisfied and the best possible solution 
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is obtained for every particular problem. If an attempt is made to 
choose suitable time-varying coefficients c(t) (kame 12) fox 
equation (3.2.3), one will find that this approach exceeds the scope 
of this thesis. Consequently, the discussion in this section will 


be limited to a stationary LHS approach. 


To start with, a numerical example will be used to reproduce 


the steps which were discussed in general in section 3.2 and 3.3. 


Suppose the following differential equation is given: 
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Hence, for the time being, two values will be taken which give a simple 
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Then (3.2.6) can be written as: 


y + 2 y + y= y + cos tor ys cos t 


The solution of the LHS is: 
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Equation (3.4.3) can be simplified by the following definition: 


Cour 
z(t) = y(t) e (Coe ey) 
so that (3.4.3) becomes: 
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If an attempt is made to solve (3.4.7) by successive 
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Although 1, (t,t) and all subsequent terms can be evaluated in closed 
form, the size of the integral makes it advisable to abandon this 


approach and to determine a bound on the modulus of 1(t,t) 
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Differentiating (3.4.9) twice with respect to t gives: 
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yields the following bound on the impulse response: 
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Considering the accuracy of the last result one might 


suggest that different numerical values for co and Cy at Sete) 


could lead to a better bound on |h(t,t)| . While the last choice 
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gave a double pole at -1 , the following decomposition, using an 
auxiliary parameter w : 


y + 2y+ Clark a y = ay Coss tata at ig +4 Cost ty 53, 4a.) 


generates two distinct poles for the LHS of (3.4.12) and might improve 
(3.4.10) for some w #0. This has been tried without success and 
(3.4.10) and (3.4.11) appear to be the best possible approximations 


that can be obtained by this - stationary LHS - approach. 


In the last example, the stationary LHS approach was used 
to determine the stability of a given system. Now the class of systems 


will be outlined which can be treated by this approach. 
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In order that (3.4.13) can be satisfied for a stationary 
LHS (3.2.3), it is necessary that (3.2.3) has poles S; in the left 
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The impulse response g(t,u) = g(t-u) , defined in (3.2.9), 


is a special solution of (3.2.3). Consequently: 
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In order that ly (t) | in (3.4.17) is bounded by a function of type 


C2. 4.13) .ei6 isprequired+that: 


eee tak Cis (3.4.18) 


Naturally, this simple derivation results in a very conser- 
vative bound (3.4.18) for the time-varying parameters. However, 
(3.4.15) represents the essential prerequisite, to apply any stationary 
LHS approach. To state the result differently, if one of the terms 


in (3.4.15) is not bounded, this approach cannot be used. 


In the remainder of this section, the stationary LHS approach 
will be applied to two further examples. In the first example, the 


following 3rd order differential equation of Popov type will be considered. 
yroytil2y+t+a(t) y=0 (ras19) 


Suppose, a range of a(t) is to be determined such that the solution 


Of (3.4.19) is stable: 
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The characteristic equation of the LHS of (3.4.20): 
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has the following solutions: Ay = Ao ON) nih 62 ys 


Consequently, the general solution of the LHS of (3.4.20) is of the 


form: 
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Now Theorem I can be applied to (3.4.26) which results in: 
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If y(t) is to be bounded by a decreasing exponential, it must be 


required that: 
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While in the last example the stationary LHS had a multiple 


root, the next example will illustrate how to use a LHS with distinct 
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Suppose the following differential equation is defined for 
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The solution of the last problem starts in the usual manner. 
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As in the case of a multiple root, it is convenient to define: 
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The integral equation (3.4.44) can be solved by differentiating twice 
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Inequality (3.4.49) is a sufficient condition that (3.4.37) is satis- 
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CHAPTER IV 


NONLINEAR SYSTEMS 


(fea Introduction 


Various aspects of linear time-varying systems have been 
dealt with in Chapter II and Chapter III. The results will now 
be applied to the original nonlinear system (1.2.1). The objective 
of this chapter is to find a sector for 9(x,t) such that the 
solution of (1.2.1) is stable in the bounded input-bounded output 


sense. 


In section 4.2 it will be shown that Theorem I can be 
applied to time-invariant nonlinear systems. A special example will 
indicate a class of systems, for which necessary and sufficient stab- 


ility conditions might be obtained. 


In section 4.3 the results from chapter III will be used 
to obtain sufficient stability conditions for a time-varying nonlinear 


system. 


In section 4.4 two integral equations will be combined in 
order to generate successive approximations for the solution of a non- 
linear time-varying system. The proof of convergence will be sketched 


briefly. 
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4.2 Nonlinear Time-invariant System 


Consider the system shown in Figure 4.A. 
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Operator 
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Figure 4.A 


The time-invariant linear operator in the forward path is 


of the following form: 
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As shown by Christensen and Trott [15], (4.2.1) is equiva- 


lent to the following integral equation: 
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Pole~shifting 


Suppose (4.2.2) is rewritten in the following form: 


E 
x(t) = f(t) + | h, (t-t) [Kexrenond (xo) cdt (4a) 
fo) 


where K is a constant. 


The unknown functions £,(t) and hy (t) can be determined in the 


following way: 
Taking the Laplace transform of. (4.2.2) 
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Therefore £,(t) and h, (t) are given by: 


a F(s) 


1+ K H(s) 
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It is interesting to note that the poles of H, (s) can 


be investigated by the customary root locus method. 


Assumption on the nonlinearity 


For later purposes it has to be required that the nonlinear- 


ity in (4.2.2) satisfies the following condition: 


6(0,t) = 0 forall ct. CED 


After these preliminaries have been dealt with, the actual 
objective of this section can be approached now. To state it in general: 
A sector for ¢(x,t) in (4.2.2) will be determined so that the solution 
of (4.2.2) is stable in the BIBO sense. The sector will be obtained in the 


following form: 
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Outline of the method 


In (4.2.5), x(t) was defined by the following integral 


equation: 
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Therefore |x(t) | satisfies the following integral inequality: 
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Similarly as in Chapter III it will be found that m1, has 
to be smaller than some positive constant 6 in order that X(t) 
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Figure 4.B 
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tnequality (4.2.12) is“1llustrated in’ Fisure 4:B°*- The 
shaded area depicts an admissible region for ec $¢(x,t) which ensures 


that the solution of (4.2.2) is bounded by at least a constant. 


In general, 6 is a function of K .. Therefore the 
admissible region for ce $(x,t) can be represented as a function 


of K , as shown in Figure 4.C, for example. 
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Example 


Suppose the following system is given: 
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Substituting the appropriate values into the following equations 
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Mo 
ik Bt 


|x(t)| <M_e CBD) 
Se 1) 


In order that |x(t)| is bounded by at least a constant, it is 


required that: 


M5 
73a 
or because of (4.2.20) 
je? - ¢ $08) <b 
which is equivalent to: 
he —ko= oC) < e ae de Gi, 2423) 


Ine Fi g@2ma.D,i nedualaty 4.2.23) is plotted with respect to 
k . For k#l1=, for example, the admissible sector for the nonlin- 


earity is shown in Fig. 4.E. 


Figure 4.D 


a, 
— 
os ; . - 7 


x ‘et at iasenos & tepsl ge x bebavod or me sal 


: z ~ om 
Sate ¢ ic 
ri: ae 
4 
(> st) & 
{> ere a3=— a 
=a 7 * 
4 
‘ ; g 3,x) 6 AN 
(£5.89) i+ *ys BeOe sng -*4 


‘ Pad 
oJ dosqao% djiw betjolq et (€8.8.8) va iisgpont,d.* .9 ff al 


<ptinon siz xo} rotose! sidleetmbe siz, .slqnsxa\102 , a ee oT 


P 


& o(xXst) | 


Figure 4.E 
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Inequality (4.2.35) is the desired result. Since K > 0 


is arbitrary, (4.2.35) can be extended to the Hurwitz sector (-1.5,~) 


4.3 Nonlinear Time-Varying System 
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First, assume that the LHS of (4.3.1) has a closed form 


solution. Then, by means of convolution, (4.3.1) can be written as: 
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Theorem I will be applied to (4.3.7) in order to get an admissible 
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Inspection of (4.3.8) shows that the last result (4.3.15) can be 
anticipated. If an attempt is made to enlarge the admissible region 
for the nonlinearity by means of pole-shifting, it is found that this 


can be done successfully for the last example. 


In general, however, the following difficulty will be 
encountered if the "pole-shifting' technique is to be carried over to 
a time-varying linear system. 
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However, if (4.3.19) satisfies the Cauchy-Lipschitz condition as 
discussed in Chapter II, then the existence and uniqueness of y(t) 
and consequently h(t,1) are guaranteed. Consequently (4.3.18) 
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Now, the approximations for |y(t)| and |h(t,t)| , 
obtained in Chapter III, can be substituted into (4.3.21). The 
remaining part of the problem can be treated in the customary manner. 
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If the terms under the double integral are bounded on [0,T] 
which is true for most physical systems, then, by Fubini's Theorem 
[16], the order of integration can be interchanged. Therefore the 


double integral becomes: 
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Equation (4.4.10) includes several special cases. 


(a) 


(b) 


(c) 


Linear time-invariant system [¢(x,t) = 0] 


For this case, (4.4.6) and (4,4.7) can be simplified to: 
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Nonlinear time-invariant system 
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(a). Consequently: 
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A sufficient condition, that {x} is a convergent 
n=0 
sequence, is: 
(x, >t) a o (x, ,t) 
| < M=cconst for all t «[0,T]. (4.4.12) 
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The proof is similar to the one given by Trott and Christensen [3] 


for time-invariant systems. To sketch it briefly: 
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Utilizing assumption (4.4.12) gives: 
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If A is chosen sufficiently large, (4.4.14) satisfies the contraction 


condition, which implies that the following series is convergent. 


z(t) = z(t) + 24(t) - a(t) + Z(t) = z,(t) eu Bane 


From™(4; 4913)¥it ‘follows’ that {x (t) } is convergent too. 
n=1 
Example 


Suppose, the first two terms of (4.4.11) are to be evaluated 


for the following system: 


x(0) 


i 
S) 


G 


x.+2 (1-.e"°) x+x+ |x| =e) (4.4.15) 


x(0) 


tt 


First, the appropriate functions y(t) and h(t,t) will 
be determined from (4.4.6) and (4.4.7) respectively. The linear 
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It should be noted that x, (t) and x, (t) depend entirely on the 
choice of (4.4.16). Therefore, there are many other possibilities to 
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CHAPTER V 


CONCLUSIONS 


The stability of time-varying nonlinear systems has been 
studied. Generally, systems are described by differential equations, 
usually not a very convenient form for stability investigations. 
Consequently, an equivalent integral equation was derived using the 


impulse response. 


The impulse response, or one-sided Green's function was 
shown to be the solution of a linear homogeneous differential equa- 
tion subject to certain initial conditions. Based on this definition, 
the existence and uniqueness of the impulse response followed from 
customary theorems of linear differential equations. Closed form 
solutions could be derived by substituting special initial conditions 
into the general solution of a linear differential equation. In 
addition, if the impulse response is to be determined on either a 
digital or an analog computer, it is simply required to solve a linear 


homogeneous differential equation subject to certain initial conditions. 


Several manipulations were applied to the integral equation 
in order to determine the stability of the appropriate solution. A 
procedure, common for every example, consisted in comparing two inte- 


gral equations, as formulated in Theorem I. 


In general, the stability condition, derived from the inte- 


gral equation, was only sufficient and not necessary. In order to 
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obtain good results, it was found advisable to avoid approximations 
when possible in the corresponding derivation. In one example, this 
rule was strictly applied which resulted in a necessary and sufficient 
stability condition. The tedious solution of this example was due 

to the fact that almost no approximation was made to simplify the 
computations. From this point of view it appears that methods which 
yield necessary and sufficient stability conditions, will remain 


limited to special classes of systems. 


As to the practical use of the methods discussed in this 


thesis, the following can be stated: 


Suppose, a particular system is completely defined and it 
is required to determine necessary and sufficient stability conditions 
only. In this case, one will likely resort to a computer, conduct a 
number of trial and error tests, and draw a stability chart in some 


parameter plane. 


However, if the problem is to design a system, the first 
preference has to be given to satisfying a certain performance criter- 
ion. According to this requirement, the control law is determined 
and the parameters are adjusted. Only then, one has to check whether 
the system exhibits 'sufficient' stability with respect to certain 
disturbances and parameter fluctuations. For this purpose, some of 


the methods, discussed in this thesis, may be adequate. 


From the integral equation, a sequence of successive approx- 


imations was generated. It could be proved that this sequence was 
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convergent under weak restrictions. Nevertheless, it should be 
pointed out here that computational difficulties often limit the 
evaluation of this series to a few initial terms. Consequently, 
this type of approximation is accurate on a 'short' time interval 


only. 
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[16] R.R. Goldberg, "Methods of Real Analysis", Blaisdell Publishing 
Company, 1964, p. 323. 
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